CERN-TH/98-150 
DTP/97/104 
May 1998 



00 

a^ 

Topology of Hadronic Flow for Higgs Production 
: at Hadron Colliders 

oo 
(N 

M. Heyssler", V. A. Khoze"''''^ and W. J. Stirling'^'' 

> 
O 

a^ 

■ " Department of Physics, University of Durham, Durham DHl 3LE, England 
lO \ 

O ■ ^ INFN - Laboratori Nazionali di Frascati, 1-00044 Frascati (Roma), Italy 
00 

0^ ■ ' TH Division, CERN, CH-1211 Geneva 23, Switzerland 

'■^i Department of Mathematical Sciences, University of Durham, Durham DHl 3LE, 

I \ England 

> 
• 1—1 

X 

M\ 

Abstract 

Hadronic radiation provides a tool to distinguish different topologies of colour 
flow in hard scattering processes. We study the structure of hadronic flow corre- 
sponding to Higgs production and decay in high-energy hadron-hadron collisions. 
In particular, the signal gg ^ H ^ bb and background gg — > bb processes are shown 
to have very different radiation patterns, and this may provide an useful additional 
method for distinguishing Higgs signal events from the QCD background. 



1 Introduction 



The distribution of soft hadrons or jets accompanying energetic final-state particles in 
hard scattering processes is governed by the underlying colour dynamics at short distances 
. The soft hadrons paint the colour portrait of the parton hard scattering, and can 
therefore act as a 'partonometer' [1T1-|12|. Since signal and background processes at hadron 
colliders can have very different colour structures (compare for example the s-channel 
colour singlet process — >• Z' — > g'g' with the colour octet process qq —>■ g* —>■ q'q'), the 
distribution of accompanying soft hadronic radiation in the events can provide a useful 
additional diagnostic tool for identifying new physics processes. 

Quite remarkably, because of the property of Local Parton Hadron Duality (see for 
example Refs. 0, 0, 0) the distribution of soft hadrons can be well described by the 
amplitudes for producing a single additional soft gluon. This distribution takes the form 
of a soft 'antenna pattern' distribution multiplying the leading-order hard scattering 
matrix element squared. Confirmation of the validity of this approach comes from recent 
studies of the production of soft hadrons and jets accompanying large Et jet and 1^+jet 



production by the CDF and DO collaborations [|T5l at the Fermilab Tevatron. 



One of the most important physics goals of the CERN LHC pp collider is the discovery 



of the Higgs boson [|T6|. Many scenarios, corresponding to different production and decay 
channels, have been investigated, see for example the studies reported in Refs. []T^ |T^. 
While final states containing leptons and photons are relatively background free, they 
generally have very small branching ratios. In contrast, the more probable decay channels 
involving (heavy) quarks have large QCD backgrounds. The question naturally arises 
whether hadronic radiation patterns could help distinguish such signals from backgrounds. 
We have in mind the following type of scenario. Suppose an invariant mass peak is 
observed in a sample of (tagged) bb events. If these correspond to Higgs production, 
then the distribution of accompanying soft radiation in the eventQ will look very different 
from that expected in background QCD production of bb pairs. One could imagine, for 
example, comparing the topologies of the hadronic flows 'on and off resonance'. 

In this study we will consider the hadronic radiation patterns for two of the standard 
Higgs processes at LHC: direct production gg —>■ H ^ bb and associated production 
qq' —>■ WH iufbb. Although the non-zero 6-quark mass is largely irrelevant when 
computing the radiation patterns, we will also consider the case when the flnal-state 
quark mass is large, so that our analysis can also be applied for example to H —>■ ti. Our 
Higgs analysis is a natural extension of the studies of Refs. where the antenna 

patterns for Z' in pp collisions and leptoquark production in ep collisions were calculated 
and shown to be different from those of the QCD backgrounds. 

The analysis presented here should be regarded as a 'flrst look' at the possibilities 



^We take this to mean the angular distribution of hadrons or 'minijets' with energies of at most a few 
GeV, well separated from the beam and final-state energetic jet directions. 
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offered by hadronic ffow patterns in searching for the Higgs. Of course, ultimately there 
is no substitute for a detailed Monte Carlo study including detector effects. However the 
results presented here indicate that the effects can be potentially large, and therefore that 
more detailed studies are definitely worthwhile. 

The paper is organised as follows. In the following section we consider direct pro- 
duction and qq decay of the Higgs boson, first for massless and then for massive quarks. 
Section 3 extends the analysis to associated production and Section 4 presents our con- 
clusions. 



2 Hadronic radiation patterns for signal and back- 
ground processes 

We begin by considering the hadronic radiation patterns for the signal {gg —^H-^ qq + g) 
and background {gg qq + g) production of a massless qq pair. The impact of non-zero 
quark masses will be considered later. The radiation pattern is defined as the ratio of the 
2^3 and 2^2 matrix elements using the soft-gluon approximation for the former. The 
dependence on the soft gluon momentum k then enters via the eikonal factors ('antennae') 

m 

PiPj 



For the QCD background process g{pi)g{p2) ^ ^'(^3)^(^4) + 9{k) we have 



1) 



9's 



\M3Wgg ^qq + g) 



with s = {pi+ p2)^ t= {Pi- Pzf 




Nr. 



Cf 
'f 



[12] + [34] 
{^[12;34i; 

^([14] + [23] -[13] -[24]) I ,(2) 



u 



{Pi - Pif, and 



[i3- kl] = 2[tj] + 2[kl] - [tk] - [tl] - [jk] - [jl] . 



(3) 



This is to be normalised by the matrix element for the leading-order scattering process 
9{Pi)9{P2) q{P3)q{PA)- 



-JM^ngg 

9; 



gg) = l(t2 + ^2) 



1 1 

Nctu 



1 1 



Cp S2 



(4) 
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The antenna pattern is then 



(5) 



\M2\'^{gg qq) 

Note that because of the non-trivial colour structure of the leading-order Feynman di- 
agrams, see Fig. |l](a), there is no simple factorisation of the eikonal factors. This is in 
contrast to the signal (Higgs) process, for which 



-2 \Mz\\gg ^qq + 

H 



n 



H 



2Nc[l2] + 2Cf[M] 



(6) 



\M2\'^{,gg qq) 

with the same momentum labeling. The two terms correspond to gluon radiation off the 
initial state gluons (colour factor and the final-state quarks (colour factor Cp)- With 
colour-singlet exchange in the s-channel (Fig. |l|(b)), there is no interference between the 
initial- and final-state emission, in contrast to the QCD background antenna pattern. 
It is this feature which will give rise to significant quantitative differences between the 
radiation patterns (see below). 
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Figure 1: The colour flow diagrams for the processes (a) gg — » qq and (b) gg ^ H ^ qq. 

The next step is to define the kinematics. The four momenta are labelled by 

a{pi) + h{p2) ^ c(p3) + d{pi) + g{k) , (7) 



where the gluon is assumed soft relative to the two large-i?T partons c and d, i.e. k <^ 
Et. Ignoring the gluon momentum in the energy-momentum constraints, working in the 
subprocess centre-of-mass frame, and using the notation p'" = {E,px,py,pz), we have 



= (-E-jnCOshr], 0, 0, i^TCOshr/) 

= {Et cosh T], 0,0, —Ft cosh 1 

= {Et cosh 7] , , Et , Et sinh 7 

P4 = {Et cosh 1], 0, —Et, —Et sinh 1]) , 



^T L-usii //, u, — J^T, — J^T 'I J ! 

;TCOsh(r7 + A?]), /c-rsin A0, fc-rcos A(/), fcTsinh(?7 + Ar])) . (8) 

This is the appropriate form for studying the angular distribution of the soft gluon jet 
relative to the large-i?T jet 3, the separation between these being parametrised by A?; 
and A0. In terms of these variables, the soft gluon phase space is 

^ ^ ^ ^ krpdkT dArj dA(j) . (9) 



(27r)3 2Ek 167r3 

We will be particularly interested in the shape of the radiation pattern as a function of the 
variables Arj and A0. Note that the direction of the soft gluon is measured with respect 
to the p3 jet. Thus for massless 2 — »• 2 scattering, coUinear singularities are located at 
Ar] = 0, A(f) = and Ar] = -21], Acj) = tt. 

We first study the QCD and Higgs radiation patterns for central qq jets, i.e. t] = 0. 
Using the kinematics of Eq. (|^) with 1] = 0, Eq. (|^) gives 

^ 4 (cosh2(Ar7) - cos2(A(/))) + Cp 

^ ^"^='^14 cosh2(Ar7)-cos2(A0) ' 

and 



2 (2 cosh^(Ar/) - cos2(A0) - 

{{4Cf - Nc) (cosh2(Ar7) - cos^{A(j)))} 



4 2Njl- cosh^ [At]) + Cf (N^ - 2) 
+ i ^ 7^4^ — (11) 



kT Nc{{4Cf - iVe) (cosh2(Ar/) - cos2(A 



Note that the radiation patterns are independent of Et- Fig. |^ shows the dependence 
of TZ^ and TZ^*^^ on At] and A0. It is straightforward to show that the patterns are 
identical close to the beam direction, 

hm 7^^.Q^^ = -iiv, , (12) 

|Ar,Hoo ki 
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gg ^ qq+g gg ^ H ^ qq+g 




Figure 2: The antenna patterns 7^«^^ of Eq. (D and = 2NJ^12] + 2C^[34] of Eq. (|) 
for the processes gg ^ qq + g and gg ^ H ^ qq + g, with r] = and = 10 GeV. The 
units of TZ are GeV~^. 

independent of A0, and close to the directions of the final-state quarks, 

lim n^'Q^D ^ I (13) 

Ar,A^~*o /cf. cosh^(Ar/) - cos2(A0) ■ 

The main difference arises from the amount of radiation between the final-state quark jets. 
To study this further we consider the distributions at the symmetric point Vc located at 
{Arj = —7] = 0,A(j) = 7r/2). This corresponds to soft gluon radiation perpendicular to 
the plane of the gg — >• qq scattering, see Fig. ||. Again using the kinematics of Eq. (||), we 
find for the QCD background process 

where the numerical value corresponds to A^^c = 3 and = 10 GeV. In contrast, for the 
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Higgs signal process we find 

7^^|,=o(Pc) = 4 (Cf + iVc) ~ 0.1733 . (15) 

There is therefore approximately 4/3 more radiation between the final-state jets for the 
Higgs production process. This is due to the absence of a colour string connecting the 
final-state quarks in the QCD background process, see Fig. 0(a). 




Figure 3: Sketch of the kinematics for the case of observing the soft gluon perpendicular 
to the event plane at the 'symmetric point' Vc = (A?7 = —77, A0 = vr/2). 

The QCD background process does, however, have colour strings connecting the 
initial- and final-state quarks, and this leads to an enhancement of soft radiation be- 
tween the jets in the plane of the scattering. We can illustrate this by considering the 
radiation patterns around the direction of the final state quark. In particular we introduce 
(as in 0) the variables (Ai?, /5), where 

Ar/ = Ai?cos/3 , A0 = Ai?sin/3. (16) 

For fixed Ai? > 0, varying 13 between and 27r describes a circle in the (A?7, A0) plane 
around the quark direction. In addition, if we fix Ai? = 7r/2 then the symmetric point 
Vc corresponds to /? = 7r/2 (or equivalently 37r/2), and the soft gluon is in the 2 — >• 2 
scattering plane for /3 = 0, vr. Figure || shows the dependence of the radiation patterns 
'j^H,QCD Qj-^ p _ as before for = final-state quarks. At /? = ti/2 ^sie 

have TZ^ > 7lP^^ , as discussed above, whereas at /5 = 0, vr we have 'RP'^^ = 7^.^.0 The 
shape of the (3 distribution therefore provides a powerful discriminator between signal and 
background. 

^In fact, the equality of the distributions at = 0, tt is true for aU AR. 
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Figure 4: The dependence of the antenna patterns TZ^ and T^Q*^^ on the angular variable 
(3 defined in Eq. (p!6D. The units of IZ are GeV~^. 

How does the interjet radiation enhancement depend on the jet rapidity rjl Again we 
consider the symmetric point located at Vc = (Ar/ = —r],A(j) = tt/2). At this point TZ^ 
is completely independent of rj, 

n"{V.) = ^{N, + CF) . (17) 

which follows immediately from Eq. (^) since [12] = [34] = at Vc- The result is 

slightly more complicated for TZ^^^. Here we find 

2 1 4C^ cosh^(r/) (N^ - 1) + iV,^ (iV, - 2Cp) 

At Vc, 7^*5^-0 jg jnaximal for r/ = with the value given in Eq. ([T^). As Ir^l ^ oo VP'~^^ 
approaches its minimum value, 

hm 7^'^^^(Pe) = ^Cf . (19) 

|r?|— ►oo Krp 

Note that in the large-77 limit the ratio IZ = IZ^ /IZ'^^^ at Vc is significantly larger than 
its value at 77 = 0: 

— , , - 7N^ + 2 

7 



7^(|^|^oo,P,) = = 3.25. (20) 

In other words, the difference in the signal and background radiation patterns at the 
symmetric interjet point increases with increasing jet rapidities. Note that the laxge-N^ 
hmits of the ratios in Eq. (|20| ) are simply 3/2 and 3, and also that 'R,{ri = 0,Vc) = 1 for 
iVc = 2. This is illustrated in Fig. ^ which shows the dependence of TZ evaluated at Vc on 
T] and Nc- 




Figure 5: The ratio TZ = TZ^ /TZ'^^^ as a function of N^. and rj at the symmetric interjet 
point At] = —1] and = 7r/2. 



2.1 Massive quarks 

So far we have only considered massless quarks. In fact for H — > bb, with <^ Mh, 
this should be an excellent approximation, since the soft gluon only 'feels' the finite b- 
quark mass very close to the jet axis, where our analysis does not in any case apply. 
Far from the jet direction, and in particular at the symmetric point Vc, the effect of the 
non-zero b mass will be negligible. The situation is however very different for the case of 
H — > tt, at Mfj > 2mt. Now mass effects are important in the radiation pattern, as we 
shall demonstrate below. 
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If we allow a finite mass for the produced quarks then the kinematics have to be 
changed accordingly. Thus we replace the kinematics of Eq. (H) by 



Pi 

P2 

Ps 

Pi 



{Eg,0,0,E,) , 

{Eg,0,pT,Eqtanhr]) , 
{Eg,0, -pT, -Egtanhr]) , 

{kr cosh(r7 + Arj), kx sin A(/), kx cos A0, kTsmh^r] + Arj)) 



[21) 



i.e. we denote the energy of the quark jets by Eg and their transverse momentum by pt- 
Thus 

Eg = cosh(r/) + pi . (22) 

We again work in the subprocess centre-of-mass frame. It is convenient to introduce the 
dimensionless variable 9 as the ratio of the final-state quark mass rrig to its energy 



e 



rrig 

Er. 



(23) 



For non-zero nig the antenna patterns receive additional contributions. For example, 
the antenna pattern of TZ^ of Eq. (R) becomes 



=7^^ -Cf[33] -Ci.[44] 



(24) 



where the massive equivalents of TZ^ and ^rP^^ are labelled with the suffix G. One effect 
of the additional terms is to cancel the final-state collinear singularities, leading instead 
to the well-known dead cone |2^ phenomenon. Using the results of Ref. we obtain 
a somewhat more complicated expression for the massive equivalent to 7l9^^ , 



-T?QCD 



= (2iV, - 2Cf + 2y) [12] + {CF-X-y) {[13] + [24]} 
+ {Cf + X -y){[U] + [23]} + 23^ [34] - C^[33] - C^[44] 



(25) 



with 



Nl 
AC. 



1 1 



U T 



1 



1 



2f^) --- -nU-,-^ +2iU-T) 



IJ2 rp2 



X 



r 1 iVci 


-1 - 


[UT Cf\ 





T^ + U'^ + 2fi- 



UT 



and 



y 



AC. 



1 

Nmr 



.UT 



Cf. 



(26) 



(27) 
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The variables T, U and /i are defined as 



T=^, f/=^, ^^=^. (28) 
P1P2 P1P2 



It is straightforward to show that the massless results are recovered in the limit mq{Q) 
0. 



2.2 Threshold behaviour (G = 1) 

We first study the behaviour of the radiation patterns Rq^^ and Rq in the threshold 
limit in which = Eg = i.e. B = 1. In fact setting 77 = we can readily 

derive the general expressions for the antennae for any value of G. Figures |^ and ^ show 
the radiation patterns for various values of B near and at threshold. Notice how the 
strong peaking structure seen in the massless case (Fig. ^ disappears as the threshold is 
approached. In fact for B = 1, the patterns do not depend on A0 at all. This can be 
seen from the analytic results. First, for B = 1 we have [34] = [33] = [44] and so, from 
Eq. (D, 

7^^=l = 2iV,[12] = ^N, , (29) 

independent of Ar] and Acj), see Fig. 0(d). 

For TZq'"^ at threshold, we first note from ( PBD that T = U = fi = ^ and thus 

From Eq. (^) we then have 

Kg™ = 2{N,-Cf + yi) [12] + (Cf - y,) {|13] + |24] + |14] + |23]} + 2(>'i - Cp)[34] 

-U'"'-^^)' '''' 

which depends on Arj but not on Acp. For |A?7| 00 T^-Q^f approaches the constant 
value 

hm 7^g2f = 2iV,[12] = ^iV, , (32) 

and becomes equal to 7^e=i' massless case. We also see from Fig. |](d) that 

T^Q^f has an absolute minimum at At] = 0, 

7^gff (Ar, = 0) = A (2iV, -C, + y,) = ^ , (33) 
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gg^qq + g = 0.7 gg^qq + g = 0.8 




gg^qq + g = 0.9 gg^qq + g = 1.0 




-2 



Figure 6: The antenna patterns IZq of Eq. (|25|) for the process gg + g with 

different values of the mass parameter G of Eq. p3|). The pseudorapidity of both quark 
jets is fixed at 77 = 0, and the transverse momentum of the soft gluon is kx = 10 GeV. In 
(d) we show the threshold result 9 = 1 {Eg = niq). The units of TZq are GeV~^. 
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gg^H^qq + g e = o.7 gg^H^qq + g e = o.8 




Figure 7: The antenna patterns TZq of Eq. (^) for the process gg H{-^ qq) + g with 
different values of the mass parameter B of Eq. p3|). The pseudorapidity of both quark 
jets is fixed at = 0, and the transverse momentum of the soft gluon is kx = 10 GeV. In 
(d) we show the threshold result = 1 {Eg = rriq). The units of TZq are GeV~^. 
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which is numerically 18% lower than the large At] value. Note the singularity in Eq. (|33| ) 
for N, = V2. 

We next consider the patterns for arbitrary rj and G. With the exception of [12] all 
antennae exhibit an r] dependence. We are again especially interested in the value of TZq 
and 7^.0*"^ at the symmetric point between the two jets at Vc = {Arj = —ri,A(f) = tt/2), 
as the massless study suggests that at this point the differences between the signal and 
background radiation patterns should be maximal. When evaluated at Vc, only [13], [14], 
[23] and [24] have an explicit 1] dependence (~ tanh(?7)), whereas 

[12] = tI, (34) 



[34] 



2-6^ 



All antennae that are rj dependent exhibit an absolute maximum at Vc of for r] —>■ 
— oo ([13], [24]) or for rj —>■ oo ([14], [23]) and vanish for rj ±oo accordingly. The fact 
that there is no rj dependence at Vc for [12], [34], [33] and [44] immediately yields (see 
Eq. (^) 

K{Vc) = ^{Nc + C^{l-e')) (35) 

for all 7], i.e. the radiation between the two jets is completely independent of their separa- 
tion in rapidity. This is illustrated in Fig. ^(a). Note that the massless result (Eq. (0)) is 
reproduced for = 0. The corresponding expression for 7^«^^ is 7] dependent and reads 

^QCD.^ . ^ 1 {4(1 - e^) + 2N^{N^ - 2)(2 - 9^)} cosh^(r/) + {2 - 9^(2 - N^)} 
® ^ N^{4CFCOsh\r])-Nc) 

(36) 

For fixed 9, TZ'q^^{Vc) always shows an absolute maximum for r/ = (see Fig. |^(b)) with 
a 9 dependence which again is maximal for the massless case 9 = 0, with the value given 
in Eq. (|T^. Once again defining the ratio of signal to background radiation patterns as 
we see that TZ^ has a local maximum at Vc, the value of which depends 
on 1] and 9, see Fig. |^(c). The value at r/ = is 

^^i^^^^ = ^)= (iV4 + 4)(92-l) + 2iV2(3-92) ' ^^^^ 

which actually shows a very weak 9 dependence. It is maximal for massless quarks 
(9 = 0) with the value (= 1.3285) already given in Eq. (pO]), and is minimal for 9 = 1 
with the value 

— N'^ -2 

ne=i{Vc,V = 0) = 4^^^^ = 1.2174 . (38) 
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Figure 8: The antenna patterns for 71q, 71q " (in units of GcV~^) and TZq at the 
symmetric interjet point Vc, for different values of the quark jet rapidity r] and the mass 
parameter ©. The soft gluon transverse momentum is taken to be kr — 10 GeV. 
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For the massless case, TZ{Vc) increased with increasing jet separation (i.e. increasing 
f]). This is again true for the massive case, as shown in Fig. |^(c). In the hmit \ri\ oo 
we find 

r T? CD ^ 4(iVc + gF(l-e^))iV,(iV,^-l) 

r,) = 2(i_e^) + Ar,2(^2_2)(2-e^) ' (39) 

which is a monotonically increasing function of O. The values at O = 0, 1 are 3.25,4.57 
respectively, for Nc = 3. 

In summary, the relative difference between the radiation patterns for the Higgs signal 
and QCD background processes is maximal at the symmetric interjet point, as depicted 
in Fig. ^. The ratio (signal/background) of the radiation patterns at this point depends 
on the rapidity of the jets and the quark mass. It is smallest (JZ = 1.33) for massless, 
central jets, and largest for massive, large-rapidity jets (JZ = 4.57). 

2.3 Radiation inside the 'dead cone' 

A final point concerns the radiation inside the dead cone of the final-state (massive) quark 
jets. In this subsection for simplicity we will only consider centrally produced jets with 
1] = — the generalisation to forward jet production is entirely straightforward. 

First we recall the result for the Higgs signal process gg —>■ H —>■ qq for massless quarks 
(see Eq. (0)), 

' cosh2(Ar7) -cos2(A0)) ' ^^^^ 

The second term is singular at the jet centre. At], A0 0, whereas the first term rep- 
resents a constant 'pedestal' of radiation from emission off the incoming gluons. In the 
massive case (G > 0), however, the singularity is removed and in fact the net contribution 
to the radiation pattern from the combination Ci?(2[34] — [33] — [44]) vanishes at the jet 
centre Vdc = {Arj = Acp = 0), hence 

7^^(Pdc,^? = 0) = -iiV,. (41) 

The corresponding result for the QCD background radiation pattern inside the dead cone 
is straightforward to calculate from the results already presented. We find, again for 

Kg°Va..-, = 0) = |^. (42) 

Interestingly, the results (|41|j4^) are independent of the quark mass, provided of course 
that rUq > 0. The effect can be seen in Figs. ^ and ^ where the value of the radiation 
patterns at their minima (i.e. inside the dead cones of the quark jets) is the same for 
all 0. The signal to background ratio in the dead cone is therefore equal to the value 
obtained at threshold and given already in Eq. 
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3 Associated Higgs Production 



Higgs production in association with a W boson qq' W* WH is a potentially 
important discovery channel at both the Tevatron and LHC colliders, especially for the 
'intermediate mass' Higgs. The non-hadronic final state WH — > £^'^77 should be rel- 
atively easy to distinguish, but unfortunately has a very small branching ratio, see for 
example the recent study in Ref. This raises the question as to whether a search in 



the decay channel 

qq' — >W* — > iiye)H{^ hi) , (43) 



might be feasible, especially with fiavour tagging of both final-state h quarks Now 
there is a potentially large irreducible background from the QCD process 

qq' > W{-^ iu^) + bb , (44) 

when Mf,i ~ Mh- The signal and background processes are illustrated in Fig. ||. 





Figure 9: Feynman graphs for the process qq' W* ih'i)H{^ qq) (associated 

Higgs production) and the background process qq' ih'i)g*{^ qq). 



We wish to study the radiation patterns for the processes (^31) and (0), in analogy 
with the gg {H —^)bb study of the previous section. We first notice that the colour 
fiows are exactly the same as those for the 2^2 scattering processes qq —>■ H ^ bb and 
qq ^ g* —>■ bb [^. We can therefore immediately write down the antenna patterns of the 
soft gluon radiation: 

7^^^ = 2C^ { [12] + [34] } - Cf [33] - Cp [44] . (45) 
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7^^^ = ^[14; 23] + 2Cp {[13] + [24]} - Cp[33] - Cp[U] , (46) 

with the momenta labelled as q{pi) + q'{p2) ^ W + h{p^) + hip^) and [14; 23] defined in 
Eq. (^). Note that the Higgs pattern is the same as for gg H bb apart from colour 
factor replacement Cp for the initial-state [12] antenna. 

In order to illustrate the quantitative differences between these radiation patterns it 
is necessary to define appropriate kinematics. Since the leading order processes are now 
effectively three-body final states, it is convenient to make some simplifying assumptions. 
Thus we assume that the H and the W are produced with zero rapidity, and that the 
b and b quarks have equal energy and have polar and azimuthal angles db and ah with 
respect to the H direction. This configuration is illustrated in Fig. ^ and corresponds to 
the four momenta 



Pi = (Vi/2, 0,0,71) 
p'i = (v^/2,0,0,-v^) 
p^H = (^H, Pth, 0,0) 

p^V = (^H^, -pT/f,0,0) 

p'^ = {Eb,PbCos{^b),pbSm{{}b)sm{ab),pbSm{{}b)cos{ab)) 

p1 = {Eb,pTH - PbCos{'db),-Pbsin{^b)sm{ab),-pbsm{{}b)cos{ab)) . (47) 

Conservation of energy and momentum gives 

EH = 2Eb= '^^}~^''' , PTH = J El - M^, Pb = JW^b. cos(^.^-^^^ 



2v^ ' ^^^^ V - , 0. V 

(48) 

The pseudorapidities and azimuthal angles of the b and b quarks are readily found to be 

^M(f>b-b) = = ^^T^i'&b-b) sm{ab-b) , (49) 

P{b,b)x 

such that ab^b = f corresponds to </)fe 5 = db,b^ 

1, f Eb+P(b;b)z \ 

Vb,b = o ■ (^°) 

2 \Eb-Pib-b)J 

The soft gluon momentum is defined relative to the 6-quark jet: 

= {kTCOsh{rib + A?]), A;rcos(06 + A0), /crsin(0f, + A0), kTsmh{rib + Ar])) . (51) 

Note that the opening angle {2i)b) between the two b quarks is a function of the partonic 
subprocess energy The dependence is illustrated in Fig. 0. Note that at threshold 
(y/§ = Mw + Mh) 2^b = 180°. 
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X 




Figure 10: The kinematics for back-to-back Higgs(-— > bb)-W production. The variables 
are defined in Eq. (|47|) . 



Let us now study the radiation patterns in more detail. We assume parameter values 
of Mh = 130 GeV, rrif, = 4.3 GeV and M^y = 80.33 GeV, and we again fix the transverse 
momentum of the soft gluon to be /ct = 10 GeV. The first thing to note is that for 
the symmetric configuration defined above, the radiation pattern for the signal process is 
independent of the azimuthal angle ab- This follows from the absence of antenna involving 
both initial- and final-state quarks in (|45|). In contrast, there is no such symmetry for 
the background process (^). 

A more striking difference is seen if we vary According to Fig. |ll] the angle 
between the final-state quarks decreases with increasing \/§ with the effect that the two 



quark jets eventually merge for large centre-of-mass energies. Figs. 12 and 13 show the 



signal and background (HB|) radiation patterns for the average value (v^ = 310 GeV) 



and for an extreme value (Vs = 14 TeV) respectively.^ The azimuthal angle ab is fixed 
at 90° in both cases. For = 310 GeV the opening angle between the b and the b 
quarks is approximately 100°. As ab = 90° the b — b plane is orthogonal to the qq' — WH 

■^Notice that at threshold, ^/Sq = Mh+Mw, the b and b are produced back-to-back, and the discussion 
is almost identical to the direct production case studied earlier, apart from colour factor differences arising 
from having incoming quarks instead of gluons. 
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Figure 11: The opening angle of the bb quark pair as a function of the partonic subprocess 
energy v^. 

scattering plane (see Fig. and thus r]b = Vl = 0- We see immediately that the main 
feature of our direct production study described earlier still holds. The most striking 
difference between the signal 71^^ and the background TZ^^ is the relative suppression 
of radiation between the 6-quark jets for the latter. There is a factor of approximately 2 
difference between signal and background radiation in the interjet region, in qualitative 
agreement with the results obtained for direct Higgs production. If we now increase the 
subprocess centre-of-mass energy the two 6-quark jets merge, forming a narrow colour 
singlet and octet state for the signal and background respectively. The situation for the 
extreme case \/l = 14 TeV is shown in Fig. |l^. Notice that for the signal process the 



soft gluon radiation becomes trapped in a very small tube. Outside the merged jets the 
radiation pattern completely flattens out. In contrast, for the background process there 
is still significant radiation between the initial- and final-state quark directions. In fact 
the distribution here is essentially identical to that for the qq' Wg process studied in 
Ref. In other words, the radiation pattern acts as a 'partonometer' in measuring 



the colour charge of the outgoing large pr partonic system. 
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4 Conclusions 



We have studied the distribution of soft hadrons or jets accompanying the production and 
qq decay of hght Higgs bosons at high-energy hadron colhders, and compared the distri- 
butions with those of the irreducible QCD backgrounds. We find significant differences 
between the signal and background distributions, which suggests that the study of the 
topology of hadron flow in such events could provide an important additional discrimina- 
tory tool. For example we have shown (see Fig. |^) that the distribution of soft hadrons 
transverse to the scattering plane in centrally produced H ^ bb events is approximately 
4/3 larger than that for QCD gg —>■ bb events with the same kinematics, while close to 
the beam axis and final-state jet directions the signal and background distributions are 
the same. The differences result from the different colour flow in the two processes. 

Although in this paper we have focused on light Higgs bosons with a dominant bb decay 
mode, we would like to make some additional remarks concerning heavier Higgs bosons. 
Consider, for example, the 'gold-plated' gg ^ H Z'^Z^ 4/^ discovery channel for a 
heavy {Mh > 2Mz) Standard Model Higgs boson. The dominant irreducible background 
comes from the qq — >■ Z^Z^ process. In the language of Section |^, the antenna patterns 
for these signal and background processes are simply 2Ai"c[12] and 2Ci;'[12] respectively. 
In other words, soft hadrons or jets with fixed transverse momentum should be uniformly 
distributed in the [t], 0) plane in both cases, but with an enhancement of 9/4 for the signal 
relative to the background. A simple on-/ off-resonance comparison should therefore show 
a significant difference. 



Our results are based on the soft-gluon/LPHD hypothesis [T^. The success of this 
approach has recently received a new quantitative confirmation from experiments at the 
Tevatron pp collider. However, it will be important to extend our work by incorporating 
a realistic Monte Carlo simulation which will allow detector effects to be included. We 
believe that the results presented in this paper make such an effort very worthwhile. 
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Figure 12: The antenna patterns (in units of GeV~^) for the signal TZq^ (Eq. (^S])) and 
the background TZq^ (Eq. (^)) for associated Higgs production at subprocess centre-of- 
mass energy y/l = 310 GeV. The directions of the incoming quarks q and q' and of the b 
and b quarks are indicated. 
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